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ABSTRACT 

We have conducted a survey of Ultra-luminous Infrared Galaxies (ULIRGs) 
with the Infrared Spectrograph on the Spitzer Space Telescope, obtaining spectra 
from 5.0—38.5 /im for 77 sources with 0.02<z <0.93. Observations of the pure 
rotational H 2 lines S(3) 9.67 /im, S(2) 12.28 /im , and S(l) 17.04 /im are used to 
derive the temperature and mass of the warm molecular gas. We detect H 2 in 77% 
of the sample, and all ULIRGs with F 60Mm > 2 Jy. The average warm molecular 
gas mass is ~2 x 10 8 M . High extinction, inferred from the 9.7 /im silicate 
absorption depth, is not observed along the line of site to the molecular gas. The 
derived H 2 mass does not depend on F 25Atm /F 60/im , which has been used to infer 
either starburst or AGN dominance. Similarly, the molecular mass does not scale 
with the 25 /im or 60 fim luminosities. In general, the H 2 emission is consistent 
with an origin in photo-dissociation regions associated with star formation. We 
detect the S(0) 28.22 fim emission line in a few ULIRGs. Including this line in the 
model fits tends to lower the temperature by ~50-100 K, resulting in a significant 
increase in the gas mass. The presence of a cooler component cannot be ruled 
out in the remainder of our sample, for which we do not detect the S(0) line. 
The measured S(7) 5.51 /tm line fluxes in six ULIRGs implies ~3x 10 6 M of 
hot (~ 1400 K) H 2 . The warm gas mass is typically less than one percent of the 
cold gas mass derived from 12 CO observations. 

Subject headings: molecular hydrogen - galaxies: ultraluminous - infrared: galax- 
ies — galaxies: AGN — galaxies: starburst 
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1. Introduction 

Ultra-luminous infrared galaxies (ULIRGs) are the most luminous objects in the local 
universe. These objects appear to be mergers of dusty, gas- rich disk galaxies (e.g., Sanders 
& Mirabel 1996; Moorwood 1996). A direct result of this dusty and large optical depth 
environment is that they radiate 90% or more of their total energy in the far-infrared, 
typically with L 8 _i oo/im > 10 12 L . The bolometric luminosities and space densities of 
ULIRGs in the local universe (i.e., z < 0.1) are similar to those of QSOs (Soifer et al. 1987; 
Sanders & Mirabel 1996). It was proposed by Sanders et al. (1988) that most ULIRGs are 
powered by dust-enshrouded QSO's in the late phases of a merger. A classic example of this is 
Mrk 231, where the detection of broad (~ 10,000 km s" 1 ) hydrogen emission lines (Arakelian 
et al. 1971; Adams & Weedman 1972) and radio jets (Lonsdale et al. 2003) provides direct 
evidence for the central active galactic nucleus (AGN). The final state of a ULIRG may be 
a large elliptical galaxy with a central massive quiescent black hole (Kormendy & Sanders 
1992). The ratio of the 25 /im to 60 fxm flux densities, as measured with the Infrared 
Astronomical Satellite (IRAS, Neugebauer et al. 1984), has been used to classify ULIRGs. 
This ratio reflects the globally averaged dust temperature of the galaxy, which is typically 
50 - 100 K. The ULIRGs are classified as either "warm" (F 2 5 Mm /F 6 o Mrrt > 0.2) and possibly 
AGN dominated, or "cold" (F 2 5 Mm /F 60Mm < 0.2), with a possible dominant contribution to 
the overall luminosity from massive star formation. Although luminous mergers are rare in 
the local universe, it has been proposed that ULIRGs may make a significant contribution 
to the star formation density at z > 1 (e.g., Elbaz & Cesarsky 2003). 

Molecular hydrogen is the primary component of the dense gas in galaxies and the most 
abundant molecule in the Universe. Stars are not only formed from molecular clouds, but the 
H 2 acts as a coolant enabling the star formation to proceed. H 2 may also act as a coolant 
facilitating the accretion of material onto a central AGN. Prior to space borne missions 
large quantities of molecular hydrogen in the central regions of ULIRGs has been inferred 
from detections of rotational transitions of 12 CO (e.g., Sanders, Scoville & Soifer 1991; and 
references for this sample in Table 4.). A much smaller, hot (~2000 K) gas component can 
be observed directly in the near- infrared through ro- vibrational H 2 lines (e.g., Van der Werf 
et al. 1993). 

In this paper, we have used observations with the Spitzer Space Telescope (Werner et 
al. 2004) to directly probe the warm (>100 K) molecular hydrogen component through pure 
rotational transitions, where the upper levels are populated via UV pumping, formation of 
H 2 in excited states, or collisional excitation. Four heating mechanisms have been proposed 
for the gas: (1) grain photoelectric heating in photo-dissociation regions (PDRs), (2) shock 
heating from outflows, SNRs and large scale streaming motions in spiral arms and bars, 
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(3) X-ray heating of gas from AGN, SNRs and cooling flows, and (4) cosmic-ray heating. 
H 2 has two independent states: ortho-H 2 , where the nuclear spins are parallel, and para-H 2 , 
where the spins are anti-parallel. The rotational states have even spatial symmetry in the 
para-H 2 and odd spatial symmetry in the ortho-H 2 . An ortho-to-para ratio of 3 is found for 
gas of typical temperature T ~ 300 K. The proton exchange mechanism is a pathway to set 
up and maintain these LTE ratios. H 2 is highly symmetric and has no permanent dipole 
moment. All the rotational transitions within the electronic ground state are quadrupolar, 
with correspondingly low spontaneous Einstein A coefficients. This makes the emission lines 
very weak, and until the advent of Spitzer, difficult to observe in all but the brightest galaxies. 
Valentijn et al. (1996) published the first detection of the H 2 S(0) 28.22 fxm transition in the 
galaxy NGC 6946, using the Short Wavelength Spectrometer (SWS) on the Infrared Space 
Observatory (ISO). Subsequently, Rigopoulou et al. (2002) published a survey of H 2 in 
12 starburst and 9 Seyfert galaxies using ISO. H 2 has been detected in only two ULIRGs, 
Arp 220 (Sturm et al. 1996) and NGC 6240 (Lutz et al. 2003). A review of the ISO 
H 2 observations is given in Habart et al. (2005). 

In order to investigate the properties for a statistically representative sample of the local 
ULIRG population, we have selected 110 ULIRGs having with redshifts between 0.02 and 
0.93 for observation with Spitzer's Infrared Spectrograph (IRS, Houck et al. 2004), as part of 
the IRS GTO program. These sources are chosen from the complete Bright Galaxy Sample 
(Soifer et al. 1987), the 1 Jy (Kim & Sanders 1998) and 2 Jy (Strauss et al. 1992) samples, 
and the FIRST/IRAS radio-far-infrared sample of Stanford et al. (2000). Armus et al. 
(2004) published the first results from this survey for Mrk 1014, Mrk 463, and UGC 5101. 
Two other papers discussing NGC 6240 and the ten ULIRGs in the IRAS bright galaxy 
sample (Soifer et al. 1987) are in press (Armus et al. 2006a & 2006b, respectively). In this 
paper we use the unprecedented sensitivity of the IRS to assess the warm H 2 component in a 
sub-sample of 77 ULIRGs. All sources were detected at 60 /zm with IRAS. Forty-eight of the 
ULIRGs in our sample are classified as "cold" (i.e., F 25Atm /F 60Mm < 0.2) and likely starburst 
dominated. The remaining twenty-nine ULIRGs are "warm" (i.e., F 25/im /F 6 o Mm > 0.2), and 
are probably powered by an obscured active nucleus. The basic properties of the sample are 
listed in Table 1. 

The Spitzer IRS wavelength coverage encompasses the pure rotational S(0) 28.22 fj,m through 
S(7) 5.51 /im transitions. The relative line strengths allow the determination of the temper- 
ature and mass of the warm H 2 gas. In §2 we detail the observations and data reduction, 
and in §3 we present the results. In general the analysis of the warm molecular gas is ap- 
proximated with a single temperature model. However, we detected the S(0) line in three 
ULIRGs, which can lower the derived warm gas temperature. Additionally, in some systems 
we detect the S(7) line, and we use this to model a hot gas component. In the discussion 
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section (§4) we investigate the warm gas mass as a function of the average global dust tem- 
perature. We also assess the warm gas fraction with respect to the cold gas mass estimates 
available in the literature. Our conclusions are presented in §5. Throughout the paper we 
adopt a flat A-dominated Universe (H G = 71 km s _1 Mpc _1 , Q M = 0.27, Q\ = 0.73, and 

= o). 



2. Observations And Data Reduction 

The data presented here were obtained using the IRS. In low resolution mode (IRS- 
LORES) there are two spectrometers with a resolving power 64 < ^ < 128. Short-low (IRS- 
SL) operates between 5.2 jum and 7.7 /im in short-low 2 (IRS-SL2) and 7.4 - 14.5 /im in short- 
low 1 (IRS-SL1). Long-low (IRS-LL) gives coverage from 14.0 - 21.3 /im in long-low 2 (IRS- 
LL2) and 19.5 - 38.0 jum in long-low 1 (IRS-LL1). The two high resolution spectrometers 
(IRS-HIRES) have a resolving power, ~ 600. Short-high (IRS-SH) encompasses the range 
9.9 - 19.6 /im and Long-high (IRS-LH) spans 18.7 - 37.2 /im. The observations were made 
in the IRS Staring Mode Astronomical Observing Template (AOT). A high- accuracy blue 
peak-up, using either a star from the 2MASS catalog (Cutri et al. 2003) or the source itself, 
was executed in order to accurately place our targets on the IRS slits. The observation log is 
given in Table 2. The IRS is fully described in Houck et al. (2004) and the observing mode 
details are presented in chapter 7 of the Spitzer Observers Manual (SOM) 1 . 

The spectral data were processed as far as the un-flatfielded two dimensional image 
using the standard IRS Sll pipeline (see the SOM). The spectra were then extracted and 
sky subtracted using the SMART analysis package (Higdon et al. 2004). To maximize the 
signal-to-noise in the final spectra the data were extracted using a column whose width 
in the cross-dispersion direction scales with the instrument point spread function. The 
spectra were flat-fielded and flux-calibrated by extracting and sky subtracting un-flatfielded 
observations of the calibration stars, HR 6348 (IRS-SL), HD 173511 (IRS-LL), and HD 
163588 (IRS-HIRES), and dividing these data by their respective templates (Cohen et al. 
2003) to generate a on-dimensional relative spectral response function (RSRF). The RSRF 
was then applied to the TARGET observations to produce the final spectra. The residual 
sky was subtracted from the IRS-LORES data using the off-source observations, which are 
part of the staring mode AOT. For the IRS-HIRES data, a spectrum of the zodiacal light 
on the date of the observation was generated using the Spitzer Planning and Observation 
Tool (SPOT. This uses the zodiacal light model of Reach et al. 2003), scaled using the IRS- 



1 http://ssc. spitzer. caltech.edu/documents/som 
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HIRES apertures, and subtracted from the spectra. The IRS-HIRES spectra were de-fringed 
in SMART. The spectra were scaled to match the 25 fim flux densities measured by IRAS. 
Sources lacking a measured IRAS 25 fim flux density were observed with the IRS' peak-up 
arrays, allowing us to scale the spectra to match the source's 22 fim flux density. The scaling 
factors, which are multiplicative, are listed in Table 2. The ULIRG IRAS 09022-36 required 
an extremely large scaling factor (~17) to bring the SH and 25 fim flux density measured 
by IRAS into agreement. We chose not to scale the data. We are confident in the detection 
of the S(l) line. However, the absolute flux calibration of this source is uncertain. 

Most of the mid-infrared emission in ULIRGs arises in a region of order one to a few 
kiloparsecs in size, which are spatially unresolved with respect to the IRS slits. 2 . To check 
this assumption we examined the IRS-SL2 data. This has the highest spatial resolution 
with 1.8" per pixel. The IRS-SL2 data was collapsed in the dispersion direction and a 
Gaussian was fit to the profile. Our calibration star, HR 6348 is fit with a Gaussian, with 
a FWHM of 2.88". Only IRAS 09022-3615, IRAS 12112+0305 and IRAS 14348-1447 show 
evidence of significant source extension i.e. FWHM > 1.5 x FWHM star . Their FWHM are 
3.0, 1.7 and 1.7 times that of the calibration star, respectively. Both IRAS 12112+0305 and 
IRAS 14348-1447 are known doubles with nuclei separated by 2.9 and 3.4 ", respectively and 
the 5" Gaussian FWHM corresponds to a linear diameter of 8 and 9 kpc, respectively. Our 
observations were centered on the NE nucleus in IRAS 12112+0305 and the SW nucleus in 
IRAS 14348-1447, as these are the dominant nuclei in the near-infrared and in CO (Evans et 
al. 2002 and Evans, Surace & Mazzarella 2000). For IRAS 09022-3615 there was a pointing 
problem and the PSF is poorly sampled. No additional correction to line fluxes has been 
made to account for the extended emission. 

Blending is a problem for the higher rotational lines, which can only be observed with 
the IRS-LORES. The S(6) 6.1 fim line, for example, can be blended with PAH emission at 
6.2 fim, S(5) 6.9 fim line can be blended with the [Ar II] 7.0 fim line and PAH emission 
at 7.0 fim, and the S(4) 8.0 fim line can be blended with PAH emission at 7.8 fim and 8.3 
fim. NGC 6240 is the only source in our sample where the emission is bright enough to be 
detected in these transitions (Armus et al. 2006a). 

The resulting S(3), S(2), and S(l) H 2 line profiles, averaged over the two slit-nod posi- 
tions, are shown in the first, second, and third columns of Figure 1, respectively. The S(0) 
detections, which indicates a somewhat cooler H 2 component, are shown in Figure 2, while 
the S(7) detections, which arise from a hotter H 2 component, are shown in Figure 3. In Table 



2 The slit sizes for the IRS modules are SL2: 3"6 x 57", SL1: 3"7 x 57", LL2: 10"5 x 168", LL1: Iff! 7 
x 168", SH: 4"7 x 11"3, LH: 11'/ 1 x 22"3 
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3 we list the line fluxes and upper- limits. We find the dominant source of uncertainty in the 
line fluxes to be the difference between the fit to the line profiles in the two slit-nod positions. 
Upper limits are calculated using the residuals (rms) from a 0-order fit to the continuum. 
We take the 3 a upper limit to be 3 x rms x FWHM. No extinction correction has been 
applied to the observed line fluxes. Determining the correct extinction factor is inherently 
difficult, and we address this in §3.2, where we discuss the location of the warm molecular 
hydrogen gas. We calculate a lower limit to the fraction of the total far-infrared luminosity, 
L IR 3 , radiated in the H 2 for those sources that have at least one measured emission line flux. 
The lower limit is the sum of the line luminosities that were fit with a Gaussian, divided by 
L IR . Typically ~ 0.01 % of the total far-infrared luminosity is radiated in the H 2 rotational 
lines. 



3. Results 

Before deriving the physical properties of the warm molecular hydrogen we will take 
a quick look at the sample to check for systematic offsets and biases. All seventy-seven 
ULIRGs were observed using IRS-LORES. However, the brighter sources were also observed 
in IRS-HIRES, fifteen with IRS-LH and thirty-seven with both IRS-SH and IRS-LH. We 
find no significant difference in the H 2 detection rate for these IRS-LORES and IRS-HIRES 
observations. For IRS-HIRES the detection rate is 35/50 (70ll2%) as compared to 51/77 
(66^9%) for IRS-LORES observations. Twenty-five sources have both IRS-LORES and IRS- 
HIRES measurement of the S(l) line. On average the IRS-HIRES measurement is 90131 
% of the IRS-LORES value. Twenty-three sources have both IRS-LORES and IRS-HIRES 
measurements of the S(3) line. Two sources, IRAS 08311-2459 and IRAS 19254-7245, have 
large line flux uncertainties. The remaining twenty-one sources have an average IRS-HIRES 
measurement that is 100^33 of the average IRS-LORES line flux. There is no evidence for 
a significant systematic difference between high- and low-resolution line flux measurements. 

The sample contains a mixture warm and cold sources. Eighteen of the twenty-nine 
warm ULIRGs were detected in H 2 (62 tl5%), whereas forty-one of the forty-eight cold 
ULIRGs were detected (85 il3%). There is thus no significant statistical difference in the 
detection rate for the two sub-populations. 

The likelihood for detecting H 2 is flux limited. There are forty sources in our sample 



3 Lir = L 8 _iooo is defined to be 4ir d^ Fir in units of L , where di, is the luminosity distance, and 
Fir = 1.8 x 10" 14 (13.48 F 12Aml + 5.16 F 25flm + 2.58 F 60flla + F wo „ m ) in W m" 2 (Sanders and Mirabel 
1996). 
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that have an IRAS 60 /xm flux density < 2 Jy. This includes the eighteen sources in the 
sample for which we do not detect any molecular hydrogen emission and twelve of the sixteen 
sources for which we only detect a single emission line. This suggests a detection limit based 
on our choice of integration time. We detect molecular hydrogen in all sources with 60 
/iin flux densities > 2 Jy, irrespective of the warm/cold classification. 



3.1. Warm H 2 Gas 

To derive the mass of warm molecular hydrogen we assume that all of our sources are 
spatially unresolved by the IRS and that the emission is optically thin. The critical densities 
of the J = 2, 3, and 4 levels are relatively low (n cr < 10 3 cm -3 ) and we assume that the 
populations are in LTE. Adopting an ortho-to-para ratio of 3, we constructed an excitation 
diagram for each source. 4 This is simply the natural logarithm of the number of molecules 
divided by the statistical weight in the upper level of each transition versus the energy level. 
If the H 2 is characterized by a single temperature, the data will lie on a straight line, with 
the excitation temperature (T cx ) being the reciprocal of the slope. The excitation diagram 
for IRAS 00188-0856 is shown in Figure 4. This is a typical example of single component fit, 
which is the case for the majority of our sources. The mass of warm H 2 can be derived from 
the line luminosity and the excitation temperature. For example, using luminosity of the 
S(l) 17.04 /iin (ortho) line, the energy of a given level is given by Ej = 85kJ(J + 1), where 
k is the Boltzmann constant. The total mass {Mxotai), is |M G , where M G is the mass of gas 
in the ortho state and M Q = m H2 N T , with m H2 being the molecular mass of H 2 and N T is 
the total number of molecules. However, N T = y^-, where Nj is the number of molecules in 

the 3 th state, Nj = ■ Here, L(J) is the line luminosity, Aj is the Einstein A coefficient, 

AEj = hz/j, where h is Planck's constant and uj is the frequency of the emission line, fj 

is the partition function for the 3 th state, fj = ^ — E , , where gj is the statistical 

£ ,/', ortho 9.7'^^ 

weight for a given state. 

We derived an excitation temperature for approximately half of our ULIRG sample (43 
sources), that had detections in at least two transitions. The mean T cx for this sub-sample 
is 336 1 15 K. Table 3 also lists an additional sixteen sources in which we only see a single 
ortho-molecular hydrogen line. We used the mean T ex to estimate the warm gas mass in 
these sources. The warm gas masses range from 10 7 to 10 9 M , with an average value of ~ 



4 The complete set of excitation diagrams is available at http://isc.astro.cornell.edu/~sjuh/ulirgs/h2/excitation- 
diagrams.html. 
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2 xlO 8 M . The warm molecular hydrogen excitation temperatures and masses are given in 
Table 4. 



3.2. Extinction and the Ortho-to-Para Ratio 

It is inherently difficult to derive the correct extinction along a given line of sight. 
Using optical depths derived in the mid-infrared alleviates the problem of extrapolating 
from a measurement made in the optical. However, the geometry of the emitting sources and 
whether the derived optical depth should be approximated as a foreground screen or material 
mixed with the warm gas, either in clumpy or uniform medium, are all highly uncertain. 
Figure 5 shows the combined low-res and high-res spectrum of Mrk 273 from 6.5 jum to 13.5 
fim. The deep silicate absorption feature corresponds to an optical depth of rg.7 MTra ~ 2, 
implying an Ay ~ 20 — 40 (Roche & Aitken 1985). For molecular hydrogen at a temperature 
of 380 K in LTE with an ortho-to-para ratio of 3, the intrinsic S(2)/S(3) line ratio is 0.56. 
If the H 2 has the same extinction implied by the depth of the silicate absorption trough, 
the line ratio will be in the range 1.05 - 1.95 when the absorbing material in a foreground 
screen, and 0.74 - 0.90 for a mixture of dust and gas. The observed S(2)/S(3) line ratio is 
0.54 1 0.08, which is consistent with the unextincted value. In Figure 6 we plot T cx versus 
the observed H 2 S(l)/S(3) line ratio, which is independent of the ortho-to-para ratio. The 
data represent the twenty-seven galaxies where at least three J < 4 lines are detected in 
the IRS-HIRES observations. NGC 6240 and IRAS 12112+0305 are not included in this 
selection, as their excitation diagrams are best described by multi-component gas, as we 
describe more fully in the next section. The figure shows that the majority of the sources 
have line ratios which are consistent with the theoretical value for no extinction. From this 
we infer that the absorbing material along the line of sight to the central AGN does not lie 
in a screen between us and the warm molecular gas. If the H 2 emission originates from a 
massive and dusty circum-nuclear torus, then the gas and dust are mixed, and the bulk of 
the emission must arise from regions with optical depths less than three, as described by 
Burton, Hollenback and Tielens (1992). Alternately, the H 2 emission may also arise in PDRs 
associated with massive star formation. Higher angular resolution observations are needed 
to measure the relative fractions emitted from a dusty torus and from PDRs. 

Using the same twenty-seven sources we can test our assumption that the ortho-to- 
para ratio is three. In Figure 7, the dashed line shows the theoretical H 2 S(l)/S(2) line 
ratio as a function of T ex , and we have over-plotted the observed line ratios. The data are 
consistent with both the assumed ortho-para-ratio of three with essentially no extinction. 
In both Figures 6 and 7, the ULIRG IRAS 15462-0450 shows the largest departure from the 
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theoretical curves, which we attribute to calibration difficulties. 

3.3. Multi- Temperature H 2 Gas 

ULIRGs are complex systems and we would not expect the H 2 gas to be at a single 
temperature. NGC 6240 (Armus et al. 2006a) is the only system in our sample where we 
detect all the S(0) through S(7) rotational lines. Its excitation diagram is shown in Figure 
8. The S(3) through S(7) lines are well fit by a hot H 2 component, with T cx = 1327 t 
36 K. The S(7) line flux implies a hot molecular hydrogen mass of (4.1 1 0.8) x 10 6 M . 
Subtracting this component's contribution from the S(0) through S(3) lines gives a second, 
cooler, H 2 component, with T ex = 292 1 6 K. The measured S(0) line flux weights the fit to 
lower temperatures and correspondingly larger H 2 masses. The warm gas mass inferred for 
the S(0) line flux is (6 1 3) x 10 8 M . 

Besides NGC 6240, we detect the S(0) line in only two other sources, IRAS 12112+0305 
and IRAS 13342+3932. Likewise, the S(7) line is detected in only six sources, NGC 6240, 
Arp 220, IRAS 06301-7954, IRAS 12032+1707, IRAS 12112+0305, and IRAS 17068+4027 
with S/N > 3. Figure 9 shows the excitation diagrams for these six sources. When the S(7) 
line is detected, a fit is first made to the S(7) and S(3) lines to determine the temperature of 
the hot H 2 component. In Arp 220, for example, the derived excitation temperature is 1435 
1 120 K. This component's contribution is then subtracted from the S(3) through S(0) lines, 
so the T cx and mass of the cooler H 2 component can be estimated. The derived T cx and 
H 2 masses are listed in Table 4. In IRAS 13342+3932 the fit is heavily weighted towards a 
low temperature, as only the S(0) and S(l) lines are detected. The resulting gas mass is an 
order of magnitude larger than the masses calculated for NGC 6240 and IRAS 12112+0305. 
A similar decrease of ~ 50-100 K in the temperature and increase in the mass by a factor of 
ten is obtained for both NGC 6240 and IRAS 12112+0305, if only the S(0) and S(l) lines 
are used in the second temperature fit. This suggests that our result for both the mass and 
temperature in IRAS 13342+3932 is strongly biased by our non-detection of the S(2) and 
S(3) lines. Using the average temperature for the sample (T cx = 336 1 15 K), we infer a mass 
of ~2xl0 8 M Q , which is probably a better estimate of the warm H 2 gas mass, and more in 
keeping with the rest of the sample. The average hot gas temperature is 1400 1 60 K, with 
an associated molecular mass of (3 1 2)xl0 6 M . The warm, T cx ~ 340 K H 2 component, 
with M H2 ~ 10 8 - 10 9 M , thus represents the bulk of the measured gas mass. 
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4. Discussion 

4.1. Warm and Cold ULIRGs 

Figure 10 shows the derived H 2 mass as a function of the IRAS 25 /im and 60 /im flux 
density ratio (F 25Atm /F 60 ^ m ). The vertical dashed line corresponds to F 25/im /F 60/im = 0.2, 
which is used to separate "cold" and "warm" ULIRGs. There are no obvious trends. That 
is, there appears to be no systematic tendency for starburst or AGN powered ULIRGs to 
have systematically larger (or smaller) masses of warm H 2 . 

Likewise, there is no clear dependence on far-infrared luminosity. Figure 11 shows the 
derived warm H 2 mass as a function of the 60 /im specific luminosity (L 60Atm ), which we 
define as 4n d^ Fgo^m in units of W Hz -1 . Here, dL is the source's luminosity distance. 
The twenty-seven ULIRGs in the plot all have z < 0.1, which ensures that the IRAS 60 /im 
flux densities correspond to nearly identical rest-frame wavelengths. A possible tendency 
for warm H 2 mass to increase with far-infrared luminosities in excess of L 60Atm = 3 x 10 25 
W Hz -1 is significantly weakened by a number of outliers with higher H 2 masses. 

The observed S(l) line intensity is typically ~10~ 4 erg cm~ 2 s _1 , which is consistent 
with the Burton, Hollenbach, and Tielens (1992) PDR models with density ~10 4 -10 5 cm -3 
and far-UV fields of ~10 3 -10 4 G . Here, G is the local interstellar far-UV radiation field, 
determined to be 1.6 x 10~ 3 ergs s _1 cm~ 2 (Habing 1968). This is consistent with PDRs 
being a large contributor to the warm H 2 emission. Of course there can also be a significant 
shock component in some ULIRGs. Constraining the contribution of shock excitation will 
require observations of ro-vibrational transitions in the near-infrared. 



4.2. Fractional Warm Gas Mass 

In Table 4 we have listed the cold molecular hydrogen mass derived from 12 CO(1-0) 
observations of thirty ULIRGs from the literature. These masses range from 10 8 to 10 11 M , 
with an average cold M H2 ~ 4 x 10 10 M Q . This is calculated using the standard Milky Way 
Mh 2 /Lco ratio of 4.6 M /K km s" 1 pc 2 determined from observations of Galactic molecular 
clouds (Solomon et al. 1987). However, as stressed by Maloney & Black (1988), because 
the CO(1-0) line intensity scales as a modest power of the gas temperature 5 , the molecular 
hydrogen masses may be overestimatedbj factors of 3-5 in regions of intense star formation or 
AGN dominance. The density and filling factor of the 12 CO(1-0) emitting regions also affect 



5 I CO oc Tfc» for T kin < 30 K, and I CO oc T kin for T ki „ > 30 K. 
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the measured CO intensity in ways that are difficult to take into account a priori. Solomon 
et al. (1997) argue that the standard conversion factor is inappropriate for ULIRGs, as 
the derived masses often exceed the dynamical masses. They argue that in the extreme 
environments found in the central regions of ULIRGs, the CO emission originates not from 
self-gravitating clouds, but from an inter-cloud medium, which is bound by the potential of 
the galaxy or molecular gas pressure, resulting in a smaller Mh 2 /Lco of ~ 1.4 M /K km s _1 
pc 2 . All of this suggests that the cold molecular masses in Table 4 should be viewed with 
considerable caution, and likely represent overestimates of the true M H2 . Indeed, use of the 
a advocated by Solomon et al. (1997) for their ULIRG sample, would both lower the cold 
H 2 mass and increase the warm H 2 gas fraction by a factor of ~3. However the warm H 2 
gas remains a small fraction of the total H 2 mass. 

We have H 2 detections for twenty-eight of these sources. Some of the sources do have 
a warm gas fraction that scales with warm gas mass, but the scatter in the data is large 
and there is no strong trend. The median fractional warm gas mass is 0.3% of the cold gas 
mass. Figure 12 shows the warm gas fraction as a function of F 2 5 Mm /F 60Atm . Rigopoulou et 
al. (2002) derived an average gas temperature of ~150 K for their sample of local starbursts 
and AGN, and found the warm gas fraction to as large as 10% in starbursts and up to 35% 
(with a large scatter) in Seyferts. In nearly half of their sample they did not detect the 
S(0) line, and the derived temperature is ~300 K. The difference between the average warm 
gas fractions in our ULIRG sample and the local starburst/AGN sample is heavily weighted 
by the adopted gas temperature, and whether or not the S(0) line is detected. It is likely 
that the warm gas fraction is similar in both. More sensitive observations of the S(0) line in 
ULIRGs are needed to fully constrain the ~150 K component. 

Figure 12 also shows that only 8/28 of ULIRGs with both cold and warm H 2 masses 
have F 25Atm /F 60Atm > 0.2, and only one has F 25Atm /F 60Atm > 0.4. As a result, we are unable to 
say conclusively if there is a tendency for the warm H 2 mass to depend on dust temperature, 
i.e., on AGN dominance. While the position of Mrk 463 in the upper-right of the figure is 
suggestive of such a trend, it may simply be an outlier, like the "cold" IRAS 12112+0305 
near the top-left. No increase in the warm gas fraction is measured as a function of either 
increasing IRAS 25 /xm or 60 /im specific luminosity. 

5. Conclusions 

We have used Spitzer's IRS to obtain spectra from 5.0—38.5 /im for a large sample of 
ULIRGs. This sample is bright (0.1 < Fqq^ < 104 Jy) and have redshifts between 0.02 and 
0.93. We have detected line emission from pure rotational transitions of molecular hydrogen 
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in 77% of the sample, and in all ULIRGs with F 60Mm > 2 Jy. The majority of these sources are 
modeled with a single temperature H 2 component, with an average excitation temperature 
(Tex) of 336 1 15 K. The ULIRGs contain large quantities of warm molecular hydrogen. 
We find an average mass 2 x 10 s M Q . Large optical depths, as inferred from the depth of 
the silicate absorption trough at 9.7 /im, are not observed along the line of site to the warm 
H 2 gas. If the gas originates in a central dusty torus, the medium must be clumpy. The warm 
H 2 mass does not appear to depend on the IRAS 25 /im to 60 /im flux density ratio, which is 
used to classify ULIRGs into "warm" (AGN dominated) and "cold" (starburst dominated) 
categories. Similarly, the warm H 2 mass does not scale with the specific luminosity at 25 /im 
or 60 /im. A possible common origin for the gas is from photo-dissociation regions associated 
with massive star formation. 

The single temperature gas model, fitted to the S(l) through S(3) lines, is a simplification 
and the gas is likely to have multiple temperature components. Detection of the S(0) line 
will typically lower the derived temperature in the single component fit by 50-100 K, leading 
to a larger estimated gas mass. This higher gas mass cannot be ruled out in the bulk of our 
sample for which we do not detect the S(0) line. A hot H 2 gas component is inferred when 
the S(7) line is detected. This line is detected in six ULIRGs, and is modeled as a T ex = 
1400 1 60 K H 2 component, with a relatively small mass (Mh 2 ~ 3 x 10 6 M Q ). 

Using cold molecular gas estimates from 12 CO observations in the literature, we find 
that the warm H 2 mass is typically less than one percent of the cold H 2 mass. The average 
warm gas fraction could be a few tens of percent, but this requires verification through the 
detection of the S(0) 28.22 /im line, which is difficult using the IRS as the line is situated 
in a low-sensitivity region of the detector array, and detection of the 12 CO line. In general, 
ULIRGs have a warm mass fraction that is comparable to that found in local starburst and 
Seyfert galaxies (Rigopoulou et al. 2002). 

This work is based [in part] on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology under 
NASA contract 1407. Support for this work was provided by NASA through Contract 
Number 1257184 issued by JPL/Caltech. 

We with to thank Paul Goldsmith, David Hollenbach, Alberto Noriego-Crespo and Terry 
Herter for useful discussions concerning the excitation of molecular hydrogen, Bill Reach for 
the Zodiacal light model, and James Houck for making the IRS possible. 
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Table 1 — Continued 
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a The luminosity distance in Mpc. 

b The 8-1000 (im far-infrared luminosity, defined as Lir, = 4-7T 
d 2 Fir in units of L , with F IR =1.8 x 10~ 14 (13.48 Fi2 Mm + 
5.16 F25/jm + 2.58 FeOfim + Fioojim), in units of W m -2 (Sanders 
and Mirabel 1996). 
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.95 










IRAS10091+47 


10:12:16.74 


46:49:42.90 


2004- 


■04- 


■19 


28. 


,5 





.87 





.88 





.84 





.85 










IRAS10398+32 


10:42:40.81 


32:31:31.00 


2004- 


■12- 


■08 


42. 


.9 





.85 





.81 





.02 





.76 











Table 2 — Continued 



Object 


RA 


Dec 


Date Obs. Exec. Time Scale Factor to IRAS 25 flux density 




(J2000) 


(J2000) 


(min) SL2 SL1 LL2 LL1 SH LH 



1RAS10565+24 


1 A 

10: 


r a 

59 


:18 


1 A 

,14 


A/i OA O/i OA 

24:32:34. oU 


AAA A 

2004- 


A r 

-05- 


1 o 

-12 


A A 

42. 


.8 


1 


A A 

.22 


1 


OA 

.32 


1 


.11 


1 


AO 

.03 


1 1 A 

1.12 


0. 


AT 

.97 


■ i > a ci i i i n i oo 


11: 


1 A 

14: 


O O 

:38 


O D 

.88 


OA /I 1 OO 1 A 

32:41:33.10 


AAA A 

2004- 


A CT 

■05- 


■11 


C A 

52. 


.0 


1 


A A 

.20 


1 


.57 


1 


A A 

,00 


1 


A 1 

,01 


1 1 A 

1.10 





.86 


■ ■ > A O 1 A A 1 O i 1 a 

1RAS12018+19 


12 


: 4: 


24. 


53 


1 A A r A O A 

19:25: 9.80 


2004- 


■05- 


■15 


51. 


.9 


1 


.06 


1 


.10 





.92 





.92 


0.87 


0. 


T/" 

.76 


TO A CM OHOO i 1 T 

IK Ab 1 20 6 2+17 


1/ 


: 5: 


-1 "7 

47. 


TO 


loiol: 8.20 


AAA A 

2004- 


A 1 

■01- 


A A 

-04 


28. 


b 


1 


A A 

.22 


1 


AT 

.2 r 


1 


AA 
.00 





A C 

.95 








TT> A C 1 A A 

!KAbl2072-04 


12 


A 

: 9: 


45. 


1 A 

.12 


E 1 1 O AA 

-5: 1:13.90 


AAA A 

2004- 


A 1 

■01- 


AC 

-06 


44. 


o 

.3 


1 


.16 


1 


AO 

.23 


1 


.11 


1 

1 


1 A 

.12 


A A A 

0.99 





.88 


■ i > * o -i oi 1 a i ao 

IHAblzllz+Uo 


12: 


13 


:46 


a tr 

.05 


A AO A 1 OA 

2:48:41.31) 


2004- 


A1 

■01- 


A A 

■04 


44. 


.0 


1 


o o 

.33 


1 


OT 

.37 


1 


.24 


1 

1 


A C 

.25 


1.12 


A 

0. 


AO 

.93 


■ ■ > a o 1 a r 1 -1 i 1 n 

1RAS125I4+10 


12 


:54 


: 0. 


.82 


1A11 1A A A 

10:11:12.40 


2005- 


■02- 


-06 


62. 


.7 





.96 





.98 





.92 





.91 


0.77 


0. 


.60 


i i ~> A O 1 O 1 A A T" A 

1RAS13120-54 


13 


:15 


: 6. 


40 


r* r* a o o o a 

-55: 9:23.30 


2004- 


■03- 


-02 


40. 


,1 


1 


.38 


1 


.49 


1 


.07 


1 


.02 


1.24 


1. 


.09 


TT7> A C1 001 O i AC 

IRAS 1.32 18+05 


1 O 

13: 


A A 

■24 


1 A 

:19 


.81 


E 07 A P C\ 

5:37: 4.60 


AAA A 

2004- 


AT 

■07- 


1 T 

■17 


51. 


1-7 

,7 





AT 

.97 





A A 

.99 





A 1 

.91 





A "1 

.91 




A 

(J. 


/"A 

.60 


TT~> A O 1 O O ,1 A i OA 

1RAS13342+39 


1 O 

13: 


36 


A A 

;24 


AT 

.07 


OA 1 T OA 1 A 

39:17:30.10 


AAA A 

2004- 


A r 

-05- 


1 A 

■12 


/"* A 

64. 


,1 


1 


AT 

.07 


1 


AT 

.07 





AT 

.97 





A /I 

.94 







.88 


IRAS 1 .3.352+04 


1 O 

13: 


36 


:51 


.15 


63:47: 4.70 


A A A C 

2005- 


AO 

■03- 


A A 

■20 


•1 T 

47. 


A 

.2 





A A 

.99 





.96 


A 




A A 

.94 


A 




A A 

.94 








TT7> A C 1 1 O ,4 C 1 i 1 A 

!KAbl3451-hl2 


1 O 

13: 


/I T 

47. 


O O 

:33 


0£? 

,36 


1 A IT A/i AA 

12:17:24. 2U 


AAA A 

2004- 


A 1 

■01- 


AT 

■07 


A A 

42. 


A 

.9 


1 


A 1 

.21 


1 


A A 

.29 


1 

1 


"1 A 

.19 


1 

1 


A 1 

.21 


1.11 


1 


A A 

.02 


TT"> AC1/iO/iO 1/1 

IKAb 14348-14 


1 A 

14: 


O T 

37 


:38 


.27 


1 £ A A A A 

-15: U:24.oU 


AAA A 

2004- 


A A 

■02- 


AT 

■07 


A A 

44. 


.4 


1 


AT 

.07 


1 


A £ 

.25 


1 

1 


1 T 

.17 


1 


1 o 

.13 


A A A 

0.99 


0. 


.81 


TT7>A01/iOTO 0/" 1 

IRAS 14378-36 


14: 


1 a 

40 


:58 


A A 

,90 


OT A OO AA 

-37: 4:33.00 


2004- 


AO 

■03- 


A A 

■02 


1 o 

43. 


2 


1 


TA 

.70 


1 


.51 


1 


A C 

.25 


1 

1 


1 A 

.19 


1.62 


1 


A A 

.02 


■ ■ > A O 1 i ro^ i 1A 

1RAS14537+19 


14 


:56 


: 4. 


43 


1 A O O A V TA 

19:38:45.70 


2004- 


-07- 


■14 


67. 


.5 





.85 





.85 





.77 





.74 




0. 


.22 


TT7> AC1 A C AO i OO 

IRAS 14548+33 


14: 


56 


:58 


A O 

.43 


OO 07 1 A A A 

33:37:10.00 


A A A C 

2005- 


AO 

■03- 


1 T 

■17 


A O 

43. 


1 





c o 

.63 





.96 


A 




A 1 

.91 





O A 

.84 








IRAS15001+14 


15: 


: 2: 


31. 


.94 


14:21:35.30 


2004- 


-06- 


■27 


56. 


.9 


1 


.09 


1 


.25 





.96 





.95 




0. 


.79 


IRAS15206+33 


15: 


22: 


:38 


.12 


33:31:36.10 


2004- 


06- 


■24 


51. 


.8 


1 


.01 


1 


.07 





.83 





.83 







.78 


IRAS15250+36 


15: 


26 


:59 


.40 


35:58:37.50 


2004- 


03- 


■04 


43. 


.4 


1 


.03 


1 


.00 





.97 





.92 


1.01 


0, 


.92 


IRAS 15462-04 


15: 


.48. 


:56 


.80 


-4:59:33.70 


2004- 


03- 


■02 


52. 


.1 


1 


.23 


1 


.35 


1 


.04 


1 


.00 


0.96 





.83 


IRAS16124+32 


16: 


14 


:22 


.11 


32:34: 3.70 


2005- 


■03- 


-19 


48. 


.0 





.41 





.92 





.88 





.87 








IRAS16334+46 


16: 


34: 


:52 


.37 


46:24:53.00 


2004- 


03- 


■04 


28. 


.8 





.83 





.35 





.65 





.63 








IRAS17068+40 


17 


: 8: 


32. 


.12 


40:23:28.20 


2004- 


■04- 


■16 


30. 


.4 





.89 





.90 





.90 





.93 








IRAS17179+54 


17: 


18: 


:54 


.23 


54:41:47.30 


2004- 


■04- 


■17 


66. 


.6 


1 


.14 


1 


.25 


1 


.08 


1 


.04 


1.18 


1 


.07 


IRAS 17208-00 


17: 


23 


:21 


.93 


-0:17: 0.40 


2004- 


■03- 


-27 


42. 


.2 


1 


.33 


1 


.35 


1 


.17 


1 


.07 


1.19 


1. 


.03 


IRAS17233+37 


17 


:25 


: 7. 


40 


37: 9:32.10 


2004- 


■06- 


■06 


67. 


.6 





.25 


1 


.12 


1 


.12 





.95 




0, 


.94 


IRAS17463+58 


17: 


:47 


: 4. 


76 


58: 5:22.10 


2004- 


04- 


■17 


32 


.8 





.96 





.96 





.91 





.92 









Table 2 — Continued 



Object 


RA 


Dec 


Date Obs. Exec. Time Scale Factor to IRAS 25 fim flux density 




(J2000) 


(J2000) 


(min) SL2 SL1 LL2 LL1 SH LH 



■ i -) * o 1 oAon i at 

1RAS18030+07 


1 o f o o r? A 

18: 5:32.50 


T O A AA 

7: 6: 9.00 


2004- 


03- 


25 


28 


5 





87 





93 





78 





75 










TT"> AOIO/i/IO i T /I 

1RAS18443+74 


18:42:54.80 


n 1 O 01 AA 

74:36:21.00 


2004- 


03- 


05 


30 


7 


1 


01 





AT 

97 





94 





95 










■ I > A 01 AOr A TO 

1RAS19254-72 


1A01 O 1 f I" 

19:31:21.55 


TO OA OO AA 

-72:39:22.00 


2005- 


05- 


30 


O T 

37 


1-7 

7 


















6 


62 


1 


31 




1n.oo.01 oc 


-o:5y:oo.3l) 


OAA A 


A C 

1)5- 


1 O 

13 


A O 


T 

7 


n 



A C 

yo 


n 



A T 

y 7 


A 

u 


O A 

o4 


A 
U 


O O 

oz 






n 



r 1 


TT> A C1 A/l CO i A A 

IRAS 19458+09 


1 A /i O 1 C TA 

19:48:15.70 


A CO CT A A 

9:52: 5.00 


OAA A 

2004- 


A -1 

04- 


1 
18 


OA 

30 


4 


1 


1 c 

15 


1 


1 A 

14 


1 


A c: 

05 


A 




96 










ttt> A conno? 1 c 

1RAS20037-15 


OA £■ 01 7fl 

20: 6:31.70 


1 CT OA O AA 

-15:39: 8.00 


onn a 

2004- 


A A 

04- 


1 O 

18 


O O 

28 


5 


1 


1 C 

16 


1 


O/" 

26 


A 




A/ 1 

96 


A 




A A 

94 










T~T> A OOAAOT AO 

1RAS20087-03 


OA 11 OO O/^ 

20:11:23.86 


O V A r A OA 

-2:59:50.80 


2004- 


05- 


14 


49 


1 





99 





96 





87 





86 









62 


■ I > A OOA1 AA A 1 

1RAS20100-41 


OA 1 O OA or" 

20:13:29.85 


A 1 /I T O /i TA 

-41:47:34.70 


2004- 


04- 


13 


44 


1 





89 





90 





86 





82 





70 





70 


1RAS20414-16 


OA H A 1 O 1 O 

20:44:18.18 


1 z" 1 /I A 1 z' 1 A A 

-16:40:16.40 


2004- 


05- 


14 


44 


1 


1 


06 


1 


07 





88 





83 









55 


■ i > A O O A r r" 1 /i c\ 

1RAS20551-42 


OA C O A/" 1 TO 

20:58:26.78 


/I O O A 1 O A 

-42:39: 1.60 


2004- 


05- 


14 


39 


7 


1 


07 


1 


11 





97 





92 









85 


TT"> A 001 0-70 i 

1RAS21272+25 


01 on on /in 

21:29:29.40 


or at crc 1A 

25:27:55.10 


OAA A 

2004- 


AC 

06- 


^ 

24 


OA 

30 


4 


1 


1 f 

16 


1 

1 


1 O 

18 


1 

1 


1 
12 


1 

1 


1 O 

13 










1RAS22491-18 


00 trl jn oc 

22:51:49.35 


1 fr r- O /) AA 

-17:52:24.00 


OAA A 

2004- 


06- 


1 

24 


44 



2 


1 


37 


1 

1 


O O 

38 


1 

1 


11 


1 

1 


11 






A 




A A 

99 


IRAS23128-59 


23:15:47.01 


-59: 3:16.90 


2004- 


05- 


11 


42 


4 


1 


59 


1 


73 


1 


43 


1 


38 






1 


11 


IRAS23129+25 


23:15:21.41 


26: 4:32.60 


2003- 


12- 


17 


40 


6 





92 





92 





92 





91 










IRAS23365+36 


23:39: 1.29 


36:21: 9.80 


2004- 


01- 


08 


43 


1 


2 


90 


2 


68 


1 


31 


1 


17 


2 


33 


1 


08 


IRAS23498+24 


23:52:26.05 


24:40:16.20 


2004- 


01- 


04 


64 


1 





97 


1 


03 





81 





84 





68 





53 


IRAS23529-21 


23:55:33.22 


-21: 2:59.60 


2003- 


12- 


16 


40 


6 


1 


00 





99 





87 





82 










MrklOW 


1:59:50.23 


0:23:40.50 


2004- 


01- 


07 


43 





1 


17 


1 


32 


1 


09 


1 


04 


1 


08 





86 


Mrk231 


12:56:14.29 


56:52:25.10 


2004- 


04- 


14 


41 





1 


20 


1 


17 


1 


03 


1 


03 


1 


18 


1 


09 


Mrk273 


13:44:42.12 


55:53:13.10 


2004- 


04- 


14 


39 


8 


1 


18 


1 


19 


1 


06 


1 


06 


1 


24 


1 


07 


Mrk463E 


13:56: 2.90 


18:22:19.00 


2004- 


01- 


07 


42 


5 


1 


00 


1 


02 





97 


1 


01 





99 





97 


NGC6240 


16:52:58.89 


2:24: 3.40 


2004- 


03- 


04 


40 








93 


1 


08 


1 


00 


1 


01 


1 


10 


1 


03 


UGC5101 


9:35:51.65 


61:21:11.30 


2004- 


03- 


23 


42 


7 


1 


11 


1 


19 


1 


13 


1 


15 


1 


19 


1 


10 



- 20 - 



- 21 - 



Table 3. Measured Line Fluxes and Limits 



Object 


S(0) 28.22 iim 
(10-2° W cm- 2 ) 


S(l) 17.04 iim 
(10- 20 W cm- 2 ) 


S(2) 12.28 (im 
(10-20 W cm- 2 ) 


S(3) 9.67 fim 
(10- 20 W cm- 2 ) 


S(7) 5.51 fim 
(lO- 20 W cm" 2 ) 


LH 2 /LiR a 

( % ) 


3C273 


< 2.2 


<3.5 


<3.2 


<3.6 


<27.3 




Arp220 


<97.0 


1.862^0.168 


0.98^0. 13 


0.73+0.02 b 


1.29^0.38 


0.006 


IRAS 00188-0856 


< 1.9 


0.043^0.005 


0.020± 0.004 


0.038±0.012 


<3.2 


0.004 


IRAS 00275-0044 


< 3.7 b 


<1.2 b 


<l.l b 


<0.5 b 


<0.8 




IRAS 00275-2859 


< 5.4 b 


<3.1 b 


<2.6 b 


<2.7 b 


<3.0 




IRAS 00397-1312 


< 1.6 


0.0291 0.004 


0.01710.002 


0.02310.006 


<8.0 


0.006 


IRAS 00406-3127 


< 5.1 b 


<1.5 b 


<1.2 b 


0.019t0.001 b 


<1.0 


0.003 


IRAS 00476-0054 




<0.1 


<0.4 b 


<0.3 b 


<0.2 




IRAS 01003-2238 


< 3.5 


0.087t0.031 


0.034±0.002 


0.057±0.007 


<2.5 


0.008 


IRAS 01199-2307 


<10.6 b 


<1.5 b 


<1.2 b 


<0.1 b 


<0.5 




IRAS 02115+0226 




<0.2 b 


<0.2 b 


0.029±0.002 b 


<0.03 


0.004 


IRAS 02433+0110 




<1.0 b 


<0.8 b 


<0.6 b 


<0.4 




IRAS 03521+0028 


< 1.3 


0.082±0.013 


0.02710.009 


0.037t0.003 


<0.3 


0.006 


IRAS 04114-5117 


< 8.7 b 


0.104t0.010 b 


<0.6 b 


0.051±0.019 b 


<0.4 


0.009 


IRAS 04313-1649 


< 6.6 b 


<0.9 b 


<0.5 b 


<0.1 b 


<0.2 




IRAS 05189-2524 


<28.2 


0.34±0.07 


0.15±0.02 


0.36^0. 13 


<34.1 


0.003 


IRAS 06035-7102 


< 6.1 


0. 417j; 0.004 


0.233^0.058 


0.34±0.10 


<11.6 


0.023 


IRAS 06206-6315 


< 2.9 


0.129^0.010 


0.050± 0.005 


0. 059^0. 019 


<1.8 


0.008 


IRAS 06301-7934 


< 9.6 b 


0.107j;0.020 b 


0.044t0.016 b 


0.072±0.017 b 


0.095±0.018 


0.022 


IRAS 06361-6217 


<14.0 b 


0.0941 0.004 b 


<2.1 b 


<0.3 b 


<2.3 


0.011 


IRAS 07598+6508 


< 3.0 


0.138t0.070 


0.037^0.006 


0.067±0.002 


<22.5 


0.011 


IRAS 08311-2459 


< 7.5 


0.388±0.030 


0.252^0.026 


0.32^0.10 


<2.6 


0.019 


IRAS 08572+3915 


<15.4 


0.11910.004 


0.05lt0.010 


0.046±0.013 


<43.0 


0.003 


IRAS 09022-3615 c 


< 3.3 


0.122±0.008 


<0.2 


<0.2 


<0.04 


0.001 


IRAS 09463+8141 


< 4.4 b 


0.108t0.055 b 


0.045t0.022 b 


0.083t0.007 b 


<0.4 


0.019 


IRAS 10091+4704 


< 4.5 b 


<0.4 b 


<0.4 b 


<0.1 b 


<0.3 




IRAS 10398+3247 




<0.1 b 


<0.001 b 


<0.004 b 


<0.2 




IRAS 10565+2448 


<10.9 


0.657±0.005 


0.256±0.009 


0.403^0.002 


<5.4 


0.013 


IRAS 11119+3257 


< 2.5 


<1.7 


<1.6 


0.054±0.007 


<9.1 


0.003 


IRAS 12018+1941 


< 2.6 


0.093^0.006 


0.030t0.003 


0.043^0.008 


<1.1 


0.010 


IRAS 12032+1707 


<10.7 b 


0.094^ 0.030 b 


<1.5 b 


0.077±0.027 b 


0.147±0.044 


0.028 


IRAS 12072-0444 


< 4.5 


0.192^0.017 


0.089^0. 010 


0.136±0.008 


<2.7 


0.020 


IRAS 12112+0305 


0.14t0.04 


0.412^0.044 


0.17310.031 


0.237t0.003 


0.41±0.08 


0.021 


IRAS 12514+1027 


< 0.7 


<0.5 


<0.5 


<0.2 


<4.0 




IRAS 13120-5453 


<30.9 


1.0471 0.034 


0.662± 0.045 


0.76310.005 


<14.6 


0.008 


IRAS 13218+0552 


< 1.8 


<1.5 


<8.9 b 


<9.4 b 


<16.9 




IRAS 13342+3932 


0. 03^0.01 


0.10210.002 


<3.9 b 


<4.2 b 


<2.6 


0.016 


IRAS 13352+6402 


< 5.2 b 


<l.l b 


<0.9 b 


0.016t0.008 b 


<1.2 


0.002 


IRAS 13451+1232 


< 4.5 


0.28310.065 


0.12910.005 


0.219±0.026 


<4.1 


0.031 


IRAS 14348-1447 


< 6.1 


0.447±0.015 


0.195±0.010 


0.240±0.021 


<3.1 


0.017 


IRAS 14378-3651 


< 6.6 


0.193t0.016 


0.07910.016 


0.12510.180 


<2.3 


0.008 


IRAS 14537+1950 




<0.1 


<0.2 b 


<0.1 b 


<0.1 




IRAS 14548+3349 




<0.2 b 


<0.2 b 


<0.1 b 


<0.2 




IRAS 15001+1433 


< 1.9 


0.077±0.007 


<2.8 b 


0.068l0.013 b 


<1.7 


0.010 


IRAS 15206+3342 


< 2.5 


0.102±0.053 b 


<4.3 b 


<2.6 b 


<1.6 


0.006 
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Table 3 — Continued 



Object 


S(0) 28.22 fan 


S(l) 17.04 /im 


S(2) 12.28 jum 


S(3) 9.67 fim 


S(7) 5.51 fim 


L H2 /LiR a 




(10~ 2 ° W cm -2 ) 


(10~ 20 W cm~ 2 ) 


(10~ 2() W cm -2 ) 


(10 -20 W cm -2 ) 


(10~ 20 W cm~ 2 ) 


( % ) 


IRAS 15250+3609 


<12.9 


0.156 + 0.014 


0.050+0.006 


0.048+0.004 


<10.4 


0.004 


TRAS 15462-0450 


< 4.1 


1 34+0 006 


037+0 001 


086+0 01 2 


<5.7 


0.010 


IRAS 16124+3241 




<0.3 b 


<0.1 b 


<0.1 b 


<0.1 




IRAS 16334+4630 


< 4.2 b 


<0.9 b 


<0.7 b 


<0.4 b 


<0.4 




IRAS 17068+4027 


<10.2 b 


0. 072^0. 018 b 


<1.7 b 


0. 052^0. 017 b 


0.087+.0.019 


0.020 


IRAS 17179+5444 


< 1.7 


0. 209+.0.007 


0.068+.0.024 


0. 14^0.02 


<1.4 


0.031 


IRAS 17208-0014 


<23.9 


0.881+0.009 


0.497+0.085 


0.57+0.11 


<8.5 


0.009 


IRAS 17233+3712 




<0.1 


<0.5 b 


<0.3 b 


<0.1 




IRAS 17463+5806 


< 3.5 b 


<0.6 b 


<0.3 b 


0. 023+0. 011 b 


<0.2 


0.007 


IRAS 18030+0705 


< 1.4 b 


0.111 + 0.008 b 


<0.8 b 


<0.4 b 


<0.2 


0.008 


IRAS 18443+7433 


<15.5 b 


<2.7 b 


<2.4 b 


0.048+ 0.007 b 


<1.7 


0.003 


IRAS 19254-7245 


< 8.5 


0.881+0.057 


0.382+0.095 


0.382+0.004 


<0.4 


0.031 


IRAS 19297-0406 


< 4.8 


0. 332+.0. 082 b 


<3.3 b 


0.196±0.008 b 


<1.7 


0.009 


IRAS 19458+0944 


<19.3 b 


0.219+_0.042 b 


<2.2 b 


0.1131 0.038 b 


<1.0 


0.009 


IRAS 20037-1547 


<11.2 b 


<5.8 b 


<4.4 b 


0.08lt0.005 b 


<4.8 


0.006 


IRAS 20087-0308 


< 2.1 


0.229+.0.014 b 


0.084+.0.033 b 


0.13j:0.04 b 


<2.2 


0.011 


IRAS 20100-4156 


< 4.6 


0.083+.0.007 


0.032+.0.001 


0.068+.0.010 


<2.9 


0.005 


IRAS 20414-1651 


< 2.4 


0.156+.0.003 b 


<1.6 b 


0.0491 0.033 b 


<0.7 


0.006 


IRAS 20551-4250 


<16.3 


0.896+.0.052 b 


<18.0 b 


0.443+.0.020 b 


<9.8 


0.013 


IRAS 21272+2514 


< 6.6 b 


<1.0 b 


<0.7 b 


0.037l0.002 b 


<0.5 


0.004 


IRAS 22491-1808 


< 7.0 


0.165t0.005 b 


<3.1 b 


0.098+.0.013 b 


<1.0 


0.006 


IRAS 23128-5919 


<12.2 


0.514t0.010 b 


<22.595 b 


0.305+.0.080 b 


<7.3 


0.009 


IRAS 23129+2548 


<10.1 b 


<1.3 b 


<l.l b 


<0.1 b 


<0.6 




IRAS 23365+3604 


< 9.1 


0.390±0.031 


0.20lt0.095 


0.31+.0.05 


<3.0 


0.015 


IRAS 23498+2423 


<0.68 


<0.5 


<0.3 


<0.2 


<2.1 




IRAS 23529-2119 




<1.4 b 


<1.0 b 


0.028l0.011 b 


<0.7 


0.010 


Mrk 231 d 


<67.5 


<31.7 


<26.8 


0.46+.0.13 


<122.3 


0.001 


Mrk 273 


<26.25 


1.0241 0.009 


0. 56^0. 07 


1. 04+.0.09 


<15.0 


0.015 


Mrk 463E 


< 7.9 


0.279+.0.037 


0.131+. 0.008 


0.28+.0.06 


<38.5 


0.017 


Mrk 1014 


< 4.2 


<2.5 


<1.5 


0.07+.0.03 


<4.6 


0.003 


NGC 6240 


0.50±0.25 


4.27lt0.344 


3.5510.26 


5.83+.0.43 


5.54+.0.50 


0.097 


UGC 5101 


<10.0 


0.496±0.047 


0.27+.0.05 


0.2810.03 


<19.5 


0.010 



a Ratio of the sum of the H2 line fluxes divided by the IR luminosity in %. The far-infrared luminosities are defined in Tabic 2. 
b Emission lines observed using IRS-LORES. 

c The absolute flux calibration for IRAS 09022-3615 is uncertain, see text 

d Both the S(l) and S(2) lines are present. However we quote upper limits as both lines are poorly calibrated. The S(l) line straddles 
orders 12 and 11 in IRS-SH and the centre of S(2) line profile is affected by bad pixel values. 
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Table 4. Properties of the Warm Molecular H 2 





Object 


T 

± ex 


Warm A/Ttt 
VvaxIIl lvlpj 2 


rv<iri Mtt 


Warm M H2 
Cold Mh 2 






I K j 


( 1 TV T \ 

(10 M Q J 


/ -i nlO a /r \ 

(10 M Q J 


( % ) 


Arp 220 


1 435 + 20 

_L tit. J '.J ZiU 


03 +0 02 

V . V • ) V7 . V7 Zj 










258 + 4 


6.6 +0.8 


3 72 Rl 


18 


IRAS 


00188-0856 


358+25 


4.11+ 1.18 


3.98 m 


0.10 


IRAS 


00397-1312 


365124 


13.161 3.69 






IRAS 


00406-3127 L 


< 336 + 15 > 


32 04+ 11 95 






IRAS 


01003-2238 


365+19 


6 68+ 2 71 






IRAS 


02115+0226 L 


< 336+15 > 


69 14+ 26 08 






IRAS 


03521+0028 


317+1 1 

OX 1 XX 


1 4 67+ 3 1 4 


5 62 m 


26 

V7 . Zj V7 


IRAS 


041 14-5 117 L 


321+26 


11.74+ 4.05 






IRAS 


05189-2524 


380+30 


2 86+1 82 


26 R2 

V7 . Zj\J 


1.1 


IRAS 


06035-7102 


363+23 


13 94+ 3 30 


4.17 ra 


0.33 


IRAS 


06206-6315 


313115 


8.141 1.76 


9.33 ra 


0.09 


IRAS 


06301-7934 L 


1329+106 
264+14 


0.3+ 0.1 
32+ 10 






IRAS 


06361-6217 L 


< 336115 > 


16.561 2.93 






IRAS 


07598+6508 


366+22 


17 53+ 9 65 


6.76 m 


0.26 


IRAS 


08311-2459 


384+22 


19 27+ 4 23 

J_ ■ / . — 1 JI . '.7 






IRAS 


08572+3915 


310118 


2.931 0.76 


0.72 ra 


0.41 


IRAS 


09022-3615 


< 336115 > 


2.611 0.48 






IRAS 


09463+8141 L 


ooo_oy 


i c i n+ o fin 
io.iu_ y.ou 






IRAS 


10565+2448 


33811 


7.051 0.10 


3.02 ra 


0.23 


IRAS 


11119+3257 


< 336115 > 


23.211 9.14 






IRAS 


12018+1941 


309111 


22.261 3.63 






IRAS 


12032+1707 L 


14711196 
309130 


0.510.5 
38120 






IRAS 


12072-0444 


34918 


19.321 2.39 






IRAS 


12112+0305 


142811049 
27115 


0.210.8 
119141 


0.49 m 


24 
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Table 4 — Continued 





Obiert 


T 

x ex 


Warm TVTtt 


Cold M H 


Warm M H2 
Cold M H2 






( If \ 


(10 IV1 J 


(10 JV1 J 


( 0/ \ 

( /o ) 


IRAS 


13120-5453 


34813 


5.291 0.24 






IRAS 


13342+3932 


141±15 
< 336115 > 


8031459 (sO) 
23.151 0.46 






IRAS 


13352+6402 L 


< 336ll5 > 


11. 10! 7.04 






IRAS 


13451+1232 


366118 


23.301 6.85 


0.51 R4 


4.57 


IRAS 


14348-1447 


32816 


19.861 1.61 


3.47 m 


0.57 


IRAS 


14378-3651 


320131 


5.971 2.42 


2.29 ra 


0.26 


IRAS 


15001+1433 L 


346140 


17.051 11.10 






IRAS 


15206+3342 L 


< 336115 > 


10.451 5.72 






IRAS 


15250+3609 


29417 


3.891 0.53 






IRAS 


15462-0450 


28515 


12.571 1.15 






IRAS 


17068+4027 L 


14161154 
292123 


0.210.2 

22l9 






IRAS 


17179+5444 


340113 


30.031 4.61 






IRAS 


17208-0014 


347113 


8.931 1.29 


3.80 m 


0.23 


IRAS 


17463+5806 L 


< 336115 > 


29.871 17.70 






IRAS 


18030+0705 L 


< 336115 > 


16.011 3.00 






IRAS 


18443+7433 L 


< 336115 > 


9.711 3.87 






IRAS 


19254-7245 


31314 


23.631 1.99 


4.57™ 


0.52 


IRAS 


19297-0406 L 


334118 


15.391 5.02 


5.25 m 


0.29 


IRAS 


i Q458-I-0944 L 


325 + 27 


1 4 Q4+ 5 80 


6 76 m 


0.22 


IRAS 


20037-1547 L 


< 336115 > 


35.921 13.50 






IRAS 


20087-0308 L 


326121 


17.491 4.69 


9.33 m 


0.1875 


IRAS 


20100-4156 


339112 


9.051 1.46 






IRAS 


20414-1651 L 


294138 


10.041 5.70 






IRAS 


20551-4250 L 


32215 


10.461 0.90 


2.82™ 


0.37 


IRAS 


21272+2514 L 


< 336115 > 


9.491 3.55 






IRAS 


22491-1808 L 


335110 


6.091 0.74 


4.27 ra 


0.14 


IRAS 


23128-5919 L 


335119 


6.041 1.36 


1.95 ra 


0.31 


IRAS 


23365+3604 


358114 


8.621 1.45 


4.68 m 


0.18 
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Table 4 — Continued 



Object 


T 

J - ex 

(K) 


Warm M Ha 
(10 7 M ) 


Cold M H2 
(10 10 M ) 


Warm Mjj 2 
Cold Mh 2 

(%) 


IRAS 23529-2119 L 


< 336115 > 


81.051 43.67 






Mrk 231 


< 336115 > 


7.971 3.73 


4.07 R1 


0.20 


Mrk 273 


37818 


6.751 0.52 


2.69 ra 


0.25 


Mrk 463E 


368121 


3.501 0.87 


0.04 R4 


8.75 


Mrk 1014 


< 336115 > 


20.381 12.20 


5.25 m 


0.39 


NGC 6240 


1327136 


0.411 0.08 








29216 


621 32 












4.37 ra 


1.43 


UGC 5101 


332111 


4.751 0.78 


0.62 ra 


0.77 



Note. — The Warm Molecular H2 Mass is derived from our data. HIRES observations are used when 
available, and the LORES observations are marked in the first column with an 'L'. The temperature is 
derived from the fit to the S(f) - S(3) data. The mass is derived using the S(I) line flux, if this line is 
not detected the S(3) line flux is used. If only a single line is detected we use the average temperature 
of 336^15 K to estimate the mass, denoted by '<>' in column 2. The standard Milky Way H2 mass- 
to-CO luminosity ratio, a =4.6 Mq /K km s _1 pc 2 has been used to estimate the cold gas mass. The 
CO references labeled RI-R4 in column 3 are: Rl: Solomon et al. 1997; R2: Rigopoulou et al. 1996; 
R3: Mirabel et al. 1990; R4: Evans ct al. 2002. However, Solomon et al. (1997) argue that a ~ 1.4 
should be used for ULIRGs. This would decrease the cold gas mass and correspondingly increase the 
warm gas fraction by a factor of ~3. The S(7) line is detected in Arp 220, IRAS 06301-7934, IRAS 
12032+1707, IRAS 12112+0305, IRAS 17068+4027 and NGC 6240, and we calculate a hot gas mass. 
The S(0) line is detected in IRAS 12112+0305, IRAS 13342+3932 and NGC 6240, which lowers the 
warm gas temperature. The S(0) line flux is used to estimate the gas mass. For IRAS 13342+3932 
we only detect the S(0) and S(l) line resulting in much lower gas temperature. We also quote the gas 
mass using the average temperature for the sample, which may be a more reasonable estimate of the 
gas mass. There is a problem with the absolute flux calibration for IRAS 09022-3615 and the mass 
could be a factor twenty higher than quoted. 
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9.56 9.619.66 9.719.76 12.20 12.24 12.28 12.32 12.35 16.93 16.98 17.03 17.08 17.13 

IRAS00188-0856 




9.61 9.64 9.66 9.69 9.71 12.20 12.24 12.2B 12.32 12.35 1 6.93 16.98 1 7.03 17.08 1 7.13 

IRAS00397-1312 IRAS00397-1312 IRAS00397- 1 31 2 




9.61 9.64 9.66 9.69 9.71 12.20 12.24 12.26 12.32 12.35 16.93 16.98 17.03 17.08 17.13 



IRAS04114-5117 



RAS04114-5117 




,.«B 




9.46 9,56 9.66 9.76 9.86 



16.83 16.93 17,03 17,13 17.23 



Fig. 1. — Pure rotational H2 emission lines for the ULIRG sample observed with Spitzer's 
IRS. The three columns show (from left to right) the H 2 S(3) 9.67 fim, S(2) 12.28 fim, 
and S(l) 17.04 /im lines after averaging over both slit positions. The data are shown as a 
histogram, while the solid line shows a single gaussian plus linear baseline fit. The vertical 
axis shows flux density in units of 10~ 18 W cm -2 /iin -1 , with rest wavelength in microns 
shown on the horizontal axis. 
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IRAS051 89-2524 



;RAS05189-2524 
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Fig. 1.— Cont. 
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IRAS09022-3615 
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Fig. 1.— Cont. 
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Fig. 1.— Cont. 
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Fig. 1.— Cont. 
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Fig. 1.— Cont. 
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Fig. 2. — Observed (histogram) and fitted (line) H 2 S(0) 28.22 jum lines detected with the 
IRS in the ULIRGs IRAS 12112+0305, IRAS 13342+3932, and NGC 6240. The spectra 
represent averages of the two slit positions. The vertical axis shows flux density in units of 
10~ 18 W cm" 2 s -1 /im -1 , while the horizontal axis shows the rest wavelength in microns. 
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Fig. 3. — Observed (histogram) and fitted (line) H 2 S(7) 5.51 /im emission lines detected with 
the IRS in the ULIRGs Arp 220, IRAS 06301-79304, IRAS 12032+1707, IRAS 12112+0305, 
IRAS 17068+4027, and NGC 6240. The spectra represent averages of both slit positions. 
The vertical axis shows flux density in units of 10~ 18 W cm -2 /iin -1 . Rest wavelength in 
microns is shown along the horizontal axis. 
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130 

1000 2000 3000 

Upper Level Energy (K) 

Fig. 4. — Molecular hydrogen excitation diagram for IRAS 00188-0856, which is typical for 
our ULIRG sample. The solid line shows a fit to the S(l), S(2), and S(3) measurements, 
which is consistent with a single temperature component at T cx = 358 1 25 K, with a mass 
of 4.1 x 10 7 M . The arrow shows the 3 a upper-limit for the S(0) transition. 
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6 8 10 12 14 

Wavelength ( Jim ) 

Fig. 5.— The combined IRS-LORES (6.5-9.8 /im) and IRS-HIRES (9.8-13.1 /im) spectrum 
for Mrk 273. The x-axis is the observed wavelength. A deep 9.7 /im silicate absorption 
trough is apparent, indicating heavy extinction (Ay ~ 20-40) towards the active nucleus. 
On the other hand, the ratio of the S(3) 9.67 /im to S(2) 12.28 /im H 2 lines infers much lower 
levels of extinction, implying that little of the absorbing material lies in a screen between 
the warm H 2 and us. 
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Fig. 6. — Molecular hydrogen S(l)/S(3) line ratios as a function of the derived temperature 
(T cx for twenty-seven ULIRGs (excluding Arp 220 and IRAS 12112+0305) with at least three 
J < 4 HIRES line measurements. The intrinsic (unreddened) ratio is shown as a dashed line. 
Extinction towards the warm H 2 component is minimal in most of the sources. 



-39- 




n L i i i i i i i i i i i i i i i i i i i i i i i i i j 

300 340 380 420 

Tex(K) 

Fig. 7. — Molecular hydrogen S(l)/S(2) line ratios as a function of the derived excitation 
temperature (T cx ) for the twenty-seven ULIRGs in Figure 8. The dashed line shows the 
expected trend for an ortho-to-para ratio of three and no extinction, which is consistent 
with the data. 
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Upper Level Energy (K) 

Fig. 8. — Molecular hydrogen excitation diagram for NGC 6240. The S(3) through S(7) 
values can be fit by with a hot (T ex = 1327 t 36 K) H 2 component. A second cooler (T ex = 
292 1 6 K) gas component is found after subtracting this component from the S(0) to S(3) 
values (re-plotted as asterisks) prior to modeling. 
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Fig. 9. — Excitation diagrams for the six ULIRGs with either S(0) and/or S(7) line de- 
tections. In five sources, the S(7) and S(3) values are fit by a single hot H 2 component, 
which is subtracted from the S(0) through S(3) values (or limits) to characterize molecular 
hydrogen gas at a cooler temperature. The derived warm (T™ arm ) and hot (T^*) excitation 
temperatures are given in the upper-right of each panel. 
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Fig. 10. — The derived warm H 2 masses for the ULIRG sample as a function of their IRAS 
25 ^m to 60 fim flux density ratios. The vertical line (F 25Mm /F 60Mm = 0.2) is thought to 
separate sources dominated by starbursts (F 25Mm /F 60/im < 0.2) from those powered primarily 
by an enshrouded active nucleus (F 25Atm /F 60Mm > 0.2). No obvious relation between warm 
H 2 mass and F 25Atm /F 60Atm is is present. 
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Fig. 11. — Warm H2 masses as a function of the IRAS derived 60 /xm specific luminosity 
(L60//m, in W Hz -1 ) for the twenty-seven ULIRGs with z < 0.1. There is no clear tendency 
for more 60 fim luminous ULIRGs to possess larger Mh 2 . 
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Fig. 12. — The ratio of the warm H2 mass to the cold H2 mass from 12 CO observations as a 
fraction of the IRAS 25 /im to 60 /zm flux density ratio. The dashed line separates "cool" 
and "warm" ULIRGs, i.e., those that appear to be powered primarily by obscured starbursts 
or active nuclei, respectively. 



